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I. INTRODUCTION 



The logical operation of a conventional digital computer does not depend on the details 
of the hardware implementation although the speed of operation and the cost of the machine 
obviously do. Conventional computers are in one particular state at any time. Furthermore 
from the point of view of programming the computer, the internal machinery of the basic 
units comprising the computer is irrelevant. This is very important because it implies that 
on a conceptual level, algorithms designed to run on a conventional computer will produce 
answers that do not depend on the hardware that is used. 

A quantum computer (QC) is very different in this respect. A QC exploits the fact that 
a quantum system can be in a superposition of states and that interference of these states 
allows exponentially many computations to be done in parallel The presence of this 

superposition is a manifestation of the internal quantum dynamics of the elementary units 
(i.e. the qubits). In other words, the quantum dynamics is essential to the operation of a 
physically realizable QC. 

The operation of an ideal QC does not depend on the intrinsic dynamics of the physical 
qubits: One imagines that the qubits are ideal two-state quantum systems that perform their 
task instantaneously and perfectly. From a theoretical point of view this situation is very 
similar to that of computers built from conventional digital circuits. However, in practice 
there is a fundamental difference. The fact that the logical operation of conventional digital 
circuits does not depend on their hardware implementation (e.g. semiconductors, relays, 
vacuum tubes, etc.) is directly related to the presence of dissipative processes that drive the 
circuits into regions of stable operation. Dissipation suppresses the effects of the internal, 
non-ideal (chaotic) dynamics of these circuits. 

The quantum dynamics of small physical devices is usually very sensitive to small per- 
turbations and this holds for qubits as well. Unfortunately, in contrast to the case of con- 
ventional circuits, dissipation usually has a devastating effect on the coherent quantum 
dynamical motion of the qubits, i.e. on the very essence of QC's. Therefore the specific 
physical realization of a QC is intimately related to the stability of its operation. 

In this paper we study the relation between the physical realization of QC's and their log- 
ical operation and explore physical mechanisms to enlarge the region(s) of stable operation. 
We demonstrate that programming a physical implementation of a QC is non-trivial, even 
if the QC consists of only two or four qubits. In most theoretical work on QC's and quan- 
tum algorithms (QA's) []5Hl^ one considers theoretically ideal (but physically unrealizable) 
QC's and therefore this problem of programming QC's (which we will call the Quantum 
Programming Problem (QPP)) is not an issue. As far as we know no experimental data has 
been published that specifically addresses this, for potential applications, very important 
and intrinsic problem of programming QC's. The aim of this paper is to investigate various 
aspects of the QPP by simulating QC hardware. 

How does a QPP reveal itself? Consider two logically independent operations (Oi and 
O2) of the machine. On a conventional computer or ideal QC, the order in which we execute 
these two mutually independent instructions does not matter: O1O2 = O2O1. However, 
it turns out that on a physically realizable QC sometimes the order does matter, even if 
there are no dependencies in these two program steps. In some cases O1O2 7^ O2O1 and 
the QC may give the wrong answers. The QPP is due to the fact that we are dealing with 
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interacting quantum mechanical objects (as communication between qubits is essential for 
computation), technical difficulties to manipulate a qubit without disturbing others and the 
fundamental physical fact that the state of a qubit cannot be frozen during the time that 
other qubits are being addressed. Also note the qualifier sometimes. There seems to be no 
general rule to decide beforehand which operation and at what stage of the QA the QPP 
leads to incorrect results. At present the only way to find out seems to be to actually carry 
out the calculations and check the results. 

It does not require a lot of imagination to see that the QPP implies that it may be very 
difficult to develop a non-trivial quantum program for a physical QC. Moreover there is no 
guarantee that a QA that works well on one QC will perform well on other QC's. Marginal 
changes in the qubit hardware may affect the interchangeability of logically independent 
operations. There are several factors that contribute to the QPP: 

1) Differences between the theoretically perfect and physically realizable one- and two- 
qubit operations, e.g. the one-qubit operations affecting other qubits and inaccuracies 
on the time-interval used to perform operations. 

2) Physical qubits cannot be kept still while others are being addressed. 

3) The effect of coupling of the qubits to other degrees of freedom (dissipation, decoher- 
ence). 

In this paper we address these issues through case studies. In Section |I| we describe the 
physical model that will be the starting point of our investigations. Our choice of physical 
models is largely inspired by NMR-QC experiments p!5|-pO|, only because other candidate 
technologies pll-p^ for building QC's are not yet developed to the point that they can 
execute computationally non-trivial QA's. 

As an example of such a QA we will take Grover's database search algorithm (see Sec- 
tion HI) and implement it on various physical models for 2- and 4-qubit QC's (Sections |3 
to |V^) . Our approach for analysing the QPP is to run Grover's QA by simulating the time- 
evolution of the physical model representing the QC. Thereby we strictly stick to the rules of 
quantum mechanics, i.e. we solve the time-dependent Schrodinger equation that describes 
the evolution of the physical apparatus representing the QC. The main vehicle for perform- 
ing these simulations is a Quantum Computer Emulator (QCE) [^]. A detailed description 
of this software tool is given elsewhere [^,^. Our work is fundamentally different from 
those of others who also address questions related to error propagation in QA's [^^-^ in 
that we execute QA's on physical models of QC hardware. The influence of non-resonant 
effects on the quantum computations using Ising spin quantum computers has been studied 



by Berman et al. |^3|. This work is similar in spirit to the one of the present paper as it 
explores the consequences of the difference between the ideal and physically realizable QC's. 
In this paper we focus on the QPP, not on methods to suppress non-resonant effects. For 



simplicity most of our calculations (Sections |^ and [V^ are done for systems at zero temper- 



ature, in the absence of interactions with the environment. Simulations of QC's operating 
at a non-zero temperature, in contact with a heat bath, are discussed in Section 0. 



3 



II. PHYSICAL MODEL 



The simplest qubit is a two-state quantum system, e.g. the spin of electrons or the polar- 
ization of photons. The basic operations in a meaningful computation are the manipulation 
of each qubit (e.g. by applying external fields) and the exchange of information between the 
qubits. In physical terms, the latter implies that there should be some interaction between 
the qubits. A non-trivial QC contains at least two qubits. It is known that the most simple 
spin-1/2 system, i.e. the Ising model, can be used for quantum computing p4|-^. In the 



presence of an external magnetic field, the Hamiltonian of the two-spin Ising model reads 

H = -J^S^SI - h{g,St + 92SI), (1) 

where J^, gi and g2 are material-specific constants, represents the applied magnetic field, 
and 3°" denotes the a-th component (a = x,y,z) of the spin-1/2 operators describing the 
nuclear spins. In this paper we use units such that h = 1. 

According to the rules of quantum mechanics, the state of the QC at time t is completely 
described by the wave function \^[t)). Executing a program on a QC is equivalent to solving 
the time-dependent Schrodinger equation (TDSE) 

z^^im) = m\m)- (2) 

For the 2-qubit QC the most general linear combination reads 



\m) = i; t)\ U) + a(T, i; t)\ Ti) + a(i, t; t)\ it) + a(t, T; t)\ TT), (3) 

where the a{bi,b2]t) are complex numbers. The probability that, upon measurement, the 
QC is in one of the four basis states | ||), . . . , | TT) is given by |a(|, |; t)p, . . . , |a(|, |; t)p 
respectively. 

The final results of a QC calculation can be read off by performing an experiment that 
measures the expectation value(s) of the spin(s). The value of a qubit is related to the 
expectation value of the 2;-component of the spin operator: 

Q,^l/2-{m\S]\m) ; j = l,2,...,iV, (4) 

where denotes the total number of qubits of the QC. In this paper we will denote the 
state of a qubit by 



|0) = IT) ; |i) = U), (5) 

and if necessary we will add a subscript, e.g. |0)i, to label the qubit. 

The Ising interaction between the spins is sufficient to implement control-NOT (CNOT) 
gates. Single qubit operations can be performed by applying additional external fields in the 
X and y direction (see below). It has been shown that any operation involving an arbitrary 
number of qubits can be written as a sequence of these elementary operations |^ , in other 
words the Ising model is a universal QC. 
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III. GROVER'S DATABASE SEARCH ALGORITHM 



Finding the needle in a haystack of N elements requires 0{N) queries on a conventional 
computer [Q. Grover has shown that a QC can find the needle using only 0{-\/N) attempts 
P.|TO[|. Assuming a uniform probability distribution for the needle, for = 4 the average 
number of queries required by a conventional algorithm is 9/4 |I8| , ^ . With Grover's QA 
the correct answer can be found in a single query [16,|18[. In general the reduction from 
0{N) to 0{\/N) is due to the intrinsic massive parallel operation of the QC. 

Experimentally this QA has been implemented on a 2-qubit NMR-QC for the case of 
a database containing four items This implementation uses elementary rotations 

about 90 degrees (clock-wise) around the x and y-axis (e.g. Xi|00) = (|00)+i|10))/y2 etc.) 
and an interaction-controlled phase shift |16 .1^. It is convenient to express these procedures 
in matrix notation, for example 



where I&1&2) = |^i)|^2) and hi = 0, 1. 
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The two qubits "communicate" with each other through the interaction-controlled phase 
shift 
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This unitary transformation can be used to implement a CNOT gate. Grover's algorithm 
can be expressed entirely by a sequence of the basic operations (6), (7), and (8). 

To see how Grover's QA works it is instructive to consider an example of a database 
with iV = 4 positions, labeled 0,...,3. Let us assume that the item we are searching for is 
located at position 2. First we put the QC in its initial state |00). Then we transform this 
state into the uniform superposition state 



|u) = -(|00) + |10) + |01) + |ll)). 



(9) 



by letting the sequence X2X2y 2-^1-^1^1 act on |00) |TI 

l^i) =X2X2F2XiXiFi|00), 

where Xj {Y j) denotes the inverse of Xj iYj). 



(10) 
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Next we apply to \u) the sequence of elementary operations |[T6|JT^j37|j3^ 

F2 = YiX,Y,Y2X2Y2l{7r), 
to encode the content of the database as 

l^2) = ^(|00) + |10)-|01) + |ll)), 



(11) 



(12) 



where we adopt the notation in which the basis states are labeled by the binary representa- 
tion of integers with the order of the bits reversed. In (12) the position of the item in the 
database (i.e. 2 in this example) is encoded by modifying \u) such that the amplitude of the 
corresponding basis state changes sign. 

The key ingredient of Grover's algorithm is an operation that determines which of the 
basis state contributes to (12) with the minus-one amplitude. In matrix notation 
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(13) 



and in terms of elementary operations 

G = XiXiY 1X2X2Y 2Y1X lY 1Y2X 2Y 2I {'n')XiXiY 1X2X2Y 2- 



(14) 



The sequence (14) is by no means unique: Various alternative expressions can be written 
down by using the algebraic properties of the X's and F's. This feature can be exploited 
to eliminate redundant elementary operations |jl8| . Starting from the uniform superposition 
state \u), one choice for the optimized sequences that implement the four different states of 
the database and Grover's search algorithm is |[l6| , p!8|j37| , |38[| 



U2 
Us 



XiY,X2Y2l{n)X,YiX2Y2l{n), 
X^Y,X2Y2Iin)X^Y^X2Y2Iin), 

XiFiX2F2/(7r)XiFiX2F2/(7r), 
XiFiX2F2/(7r)XiFiX2F2/(7r), 



(15a) 
(15b) 
(15c) 
(15d) 



where the [/„ correspond to the case where the needle is in position n. 

The sequences of unitary transformations, e.g. Xi,l2,etc., are easily emulated on an 
ordinary computer. We use our software package, Quantum Computer Emulator (QCE) 
|37| , to perform these calculations. An example will be shown in Appendix 0. On an ideal 
QC, sequences (0) return the correct answer, i.e. the position of the searched-for item. 
This is easily verified on the QCE by selecting the elementary operations (called micro 
instructions on the QCE) that implement an ideal QC. 



IV. TWO-QUBIT QC'S 

The energy-level structure of the nuclear spins of molecules such as deuterated cytosin 
||T5| , p^ , ^ and carbon-13 labeled chloroform |[T7|JT8[] can be described by model (1) and hence 
they can be used as physical realizations of 2-qubit QC's. 
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A. Resonant pulses 



NMR-QC experiments on carbon-13 labeled chloroform use resonant pulses to ma- 



nipulate the quantum state of the nuclear spins of the and ^'^C atoms. In the presence of 
a static magnetic field along the +z direction this NMR-QC system is described by (1) with 



h^gi/2'K ^ 500MHz, Kg2l1i[ ^ 125MHz, and J^jl-n -215Hz A detailed account 

of simulations for this case have been published elsewhere [^,^]. Simulations for physical 
model (1) confirm that sequences (|T5]) yield the correct answers for the database search 
problem |3^. However, we also demonstrated that the outcome of these calculations may 



be very sensitive to the order in which logically independent operations are carried out [BS .. 
Some of these results are reproduced in Table | (first six rows, for details see []37| , |38| ). 

The results marked with a tilde are obtained by using a logically identical but physically 
different uniform superposition, i.e. 



= AiAiriXaXsralOO). (16) 

On an ideal QC, = \u) but on a physical realizable machine this is unlikely to be the 
case. In an experiment it is simply impossible to freeze spin 2 (1) during the time that 
resonant pulses are being applied to spin 1 (2). Unless the length of these pulses is chosen 
judiciously, the wave function will acquire an additional phase. The corresponding unitary 
transformation does not necessarily commute with the operations that follow, potentially 
leading to an incorrect final result (as shown in Table |), as we pointed out in a previous 
paper 
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For the two-spin system (16) one may optimize the pulse durations such that the effect 
of these phase errors yields qualitatively correct answers. A basic step is to make the pulse 
lengths commensurate with all relevant time scales E3,H6|. The main idea of the method 



suggested in |^ and demonstrated in is to choose the parameters of the electromagnetic 
pulse such that the non-resonant spin still rotates but returns to its initial position at the 
end of the pulse . The results of the calculations are shown in Fig. ^ and summarized in 
Table | (primed symbols) [Q. In the numerical calculations we used J = — 10~^, h^gi = 1, 
hzg2 = 0.25 , {hyQi = 0.025, hyg2 = 0.0625, ^p^igg = 47r), for the pulses on the first spin 
and {hyQi = 0.05, hyg2 = 0.0125, ^p^jgg = Svr) for the pulses on the second spin, all values 
in dimensionless units. It is clear that optimization has the desired effect on the sequences 
that operate on \u') (symbols with a tilde). 

In spite of this optimization. Table | shows that there are significant quantitative differ- 
ences between the theoretically exact results (rows (1,2)) and those obtained by simulating 
a physical model of a QC (e.g. rows 7 to 10). Even if the pulse length is taken to be com- 
mensurate with the relevant time scales of the QC, changing the state of the qubits by way 
of resonant pulses yields quantum states that are different from those obtained by means 
of the unitary transformations used in the analysis of the ideal QC. This is because spin 2 
also interacts with the field applied to spin 1 and vice versa. With each program step, the 
non-ideal unitary operation may or may not result in the proliferation of errors. These errors 
are systematic (there is no "random" error source in our calculations) and directly linked 
to the structure of the QA. This is a clear case of a QPP, although we managed to let the 
different QA's produce the correct answer. Note that the QPP cannot be solved by means 
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of error correction [^]: The operations on the extra qubits required for error correction will 
suffer from exactly the same QPP. 



B. Stability of Grover's quantum algorithm 

So far we studied the stability of quantum algorithms by perturbing the input to the 
database encoding part of the algorithm. In this subsection we will study the QPP of the 
database query part of Grover's search algorithm (the operation G, see (14)) itself. First we 
will assume that the input provided to G is exact (i.e. of the form (10) for example) and we 
will compare the output of logically identical but physically different implementations of G. 
As examples we take the original sequence 

G = XiXiYiX2X2Y2YiXiYiY2X2Y2l{n)XiXiYiX2X2Y2, (17) 

and two, logically identical, sequences 

G = XiX,YiX2X2Y2YiXiY2X2YiY2l{n)XiXiYiX2X2Y2, (18) 



G = XiXiY,X2X2Y2Y2X2Y2Y^X,YJ{7i)X2X2Y2X^X^Y^. (19) 

Note that on purpose we did not "optimize" these sequences by using e.g. XiXi = 1. On 
an ideal QC we have 



G\^2) = G^\^2) = G\^2) = G=^|^2) = G\^2) = &\^2) = |01), (20) 

providing another test of the stability of the query operation G on a physical QC. Table |l| 
contains the numerical results obtained by running the sequences (17), (18), and (19) on 
the QCE using the exact state |\E'2) (see (14)) as input. In the case of the NMR-like 
QC optimized resonant pulses were used. The ideal QC performs as expected but the 
physical implementation (symbols with a prime) does not. In fact even one application of 
G apparently returns an answer that is close to being useless {Q[ ~ 0.5). As the three 
sequences (17), (18), and (19) are logically identical this is a clear case of a QPP. 

The occurrence of a QPP seems to be a generic feature of QA's running on QC's. There- 
fore it is of interest to try to quantify the QPP. We now describe a simple procedure for this 
purpose, using G and the case where the item is located in position (the exact input state 
being |\E'o)) as an example. In general there are two sources of errors in this calculation: The 
input |\E'') to G and G itself, the latter depending on the particular hardware implementation 
of the QC. As before we will use the resonant pulse technique in our numerical experiments. 

We write |\E'') as 

1^') = aol^o) + ail^i) + 021^2) + asl^s), (21) 

where the amplitude ao can always be taken real (— 1 < ao < 1) and is chosen at random. 
The other three complex coeffients are chosen randomly too, subject to the constraint + 
|tt2p + Ic^sP = 1 Ic^oP- The real variable x = (^'|^o) = c^o is a measure for how much 
the input state deviates from the exact reference input I^E'o)- On an ideal QC, Gj^E^o) = 
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1 00). Thus we can use the state |00) as reference to determine how much the output state 
|0) = G'|\E'') deviates from the exact answer. We quantify this deviation by the variable 
y = \mG^')\. 

The result of a numerical experiment using 20000 random input states |\E'o) is shown in 
Fig. |1|. Plots for the three other cases are nearly identical. We classify input-output pairs 
as "good" or "bad" as follows. First we choose a confidence level 0<c<l(c = 0.7 for the 
data shown in Fig. |I|). A particular input-output pair is considered to be good if a;^ > c and 
y"^ > c. In Fig. |l] the good (bad) pairs are shown by black (gray) markers. 

At a fairly low confidence level of c = 0.7, the region of stable operation of the G operation 
is rather small. This corroborates our earlier finding that successive applications of G, e.g. 
G^, rapidly drive the system into a region of instability. In general, quantum systems are 
very sensitive to noise and become more sensitive as the number of operations on qubits is 
increased @]. In the absence of dissipation, it is easy for the system to leave the relatively 
small manifold of good input states, a characteristic feature of almost chaotic dynamics. 



C. Hard nonselective pulses 

Another physical implementation of a 2-qubit QC employs the nuclear spins of two 
spin-1/2 nuclei in a solution of cytosine in D2O JI3tp!B| , P?[ . This system can also be described 
by Hamiltonian (1). In the NMR experiments hard nonselective pulses are used to address 
the qubits. In this section we study the stability of QC operation for this physical realization 
of a QC. 

As usual it is expedient to transform to another frame of reference that rotates with a 
constant frequency. This is accomplished by substituting in the TDSE 

|$(t)) = e**^-(si+<;2)(5f+5|)/2|^(^^^)^^ ^22) 

The time evolution of {"^{t)) is then governed by the Hamiltonian 

H = -J^SISI - ^Sl + l^l, (23) 

where VL = hz{gi —92)- Guided by experiment ||15| , ^ in our numerical work we will set 
VL/2 = 1 and J^/vrfi = —0.01887 (in dimensionless units). 

We now consider the time evolution of the two spins when we apply a static magnetic 
field hj. along the x-axis. The Hamiltonian in the laboratory frame is given by 

H = -J.SISI - KigiSt + 92S^) - h^igiSl + g^S^^), (24) 

and the corresponding expression in the rotating frame of reference reads 

- K{9iSl + 928^) cosnt - K{9,Sf + 92SI) sinOt. (25) 

If the duration of the pulse is short, i.e. VLt <^ 1, it is a good approximation to drop the 
time-dependence in (25) and we obtain 

H = + - K{9iS'^ + 928^), (26) 
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where we used the fact that since | <^ |f^/2|, for short pulses the effect of the spin-spin 
interaction on the time evolution is small and may be neglected. It is instructive to compute 
the time evolution of the spins, initially in state |0) (|0) = | |) by convention), under these 
circumstances. A straightforward calculation yields 



where Xj = + hlgj)^/^. Expression (27) shows that for hard pulses Hh^gjl > ^/2) 

and Xjt = 71, the effect of the pulse is to change qubit j from |0) to approximately |1). Note 
however that a sequence of such pulses can never turn |0) into |1) exactly and that both 
spins are affected by the pulse. 

In the rotating frame of reference the two spins rotate around the z-axis in the opposite 
sense. This fact can be exploited to perform operations that leave one spin intact while 
flipping the other spin P9| , p2| , p3[ . For instance, to rotate spin 1 around 90 degrees (clockwise) 
about the x-axis and without changing the state of spin 2, an operation we call Xi, we can 
use the sequence of pulses 

X^ = Xi,Y,^,Zl,Y,^„ (28) 

where X[ 2 denotes a 45-degree pulse around the x-axis acting simultaneously on both spins 
(see (26)), Yi^2 a 90-degree pulse around the y-axis, and Z[ 2 represents a pulse during which 
the spins evolve in time according to (26) and make a rotation of 45 degrees around the 
2;-axis, (clockwise in the case of spin 1, anti-clockwise in the case of spin 2). There are many 
sequences of two-spin pulses that yield Xi (or, as a matter of fact, any other rotation). For 
example the inverse of X^, X^, can be obtained directly from (28), i.e. 

^1=^1,2^1,2^1,2-^1,2) (29) 

or can also be written as 

^1 = -^1,2^1,2^1,2^1,2- (30) 

For our numerical calculations we have chosen to work with the sequences 

X-^ = -^1,2^1,2^1,2^1,25 (31a) 

^1 = -^1,2^1,2^1,2^1,2' (31b) 

X2 = -^1,2^1,2^1,2^1,2) (31c) 

-^2 = -^1,2^1,2^1,2^1,2) (31d) 

^1 = ^1,2-^1,2^1,2-^1,2) (31e) 

^1 = ^1,2-^1,2^1,2-^1,2) (31f) 
^2 = ^1,2-^1,2^1,2-^1,2) (31g) 
^2 = ^1,2-^1,2^1,2-^1,2) (31h) 

to perform the single-qubit operations that are used in the implementation of Grover's 
search algorithm. As in the case of resonant pulses the natural time-evolution of the system 
provides the tool to carry out the conditional phase-shift operation /(vr). 
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The numerical values for the qubits 1 and 2 in the final state of the machine, obtained 
by running the QCE for this particular hardware realization of a QC and QA are are given 
in Table Q 

From a comparison with the data of Table | we conclude that this hardware realization 
of a 2-qubit QC seems to be much less sensitive to QPP errors. As a matter of fact, in 
our numerical experiments (data not shown but included in the QCE download) we found 
none. The main reason for the good performance of this QC can be traced back to the fact 
that the QC contains two spins with exactly the same precession frequency, precessing in 
opposite direction. It is clear that it will become rather cumbersome if not impossible to 
design a A^-qubit QC (if N 2) using this approach: We would have to solve a complicated 
optimization problem that for each QA we would like to run on the QC. 



D. Reducing systematic errors 

For QC hardware to be of practical use, a basic requirement is that logically identical 
but physically different implementations of the same QA yield the same answers. The 
examples discussed above show that optimization of the pulses is essential to achieve this 
goal. For 2-qubit QC's, the optimization problem is relatively easy to solve because only two 
frequencies are involved. Alternatively one may choose to use the same duration for all pulses 
and optimize the strengths of the pulses. However, in the case of A^ ^ 2 qubits, choosing 
pulse durations commensurate with A^ different frequencies will considerably slow down the 
operation of the whole machine. In addition, the longer the pulses, the more important the 
effect of imperfections of the pulses becomes. Although there are ingenious techniques to 
optimize pulses in this respect [^^4| -|56[|, it may well be that in order to successfully run 



a QA on A^-qubit QC hardware one first has to solve a rather complicated optimization 
problem and then simulate the QC by running the QA on a conventional machine. 



V. CAN DISSIPATION REDUCE THE QPP ? 

Elsewhere we reported on the effect of dissipation on the quantum dynamics of nanoscale 
magnets ||57| , |58[| . In these systems dissipation usually causes decoherence. Decoherence 



limits the time for performing quantum computations. Above we have demonstrated that 
programming non-ideal QC's is difficult due to intrinsic instabilities of the physical device. 
Therefore we wonder if dissipation can provide the stabilizing processes, perhaps at the cost 
of reducing the time of coherent quantum operation. 

Let us consider the procedure to generate the uniform superposition state from the state 
with all spins up. Above we made use of sequences such as XXY but it is easy to see that 
from a programming-point-of-view Y is just as good: In a classical picture both sequences 
change the direction of a spin from the z to the x-axis. Thus, on an ideal QC we may 
omit XX if we want. However, in an experiment (in which there always is some amount 
of dissipation), XX can help to stabilize the direction of the spin once it points in the x- 
direction. Although in the pure quantum case XX only causes the spin to precess around 
the X-axis, the damping of this motion by dissipation will let the spin relax to the x-direction. 
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equation |^| 



We have studied the effects of dissipation on quantum computations using the master 
^P^^^ -i[H, p{t)] - A ([C, Rp{t)] + [C, i?p(t)]t) , (32) 



dt 

that describes the time-evolution of the reduced (to the subspace of the qubits) density 
matrix p{t). The operator R is defined as 

{k\R\m) = CiEk - Era)Np{Ek - E„,){k\C\m), (33) 

where H\k) = Ek\k), H being the Hamihonian describing the QC. We take for the spectral 
density of the bosonic thermal bath C{E) = IoE'^sign{E) , i.e. an super-Ohmic reservoir, and 
Np{E) = {e^^ — 1)~^ denotes the boson occupation number [Q. The operator C specifies 
the coupling between the spins and the bosons, e.g. C = J2j{Sj + Sj)/2,. We set the inverse 



temperature P = 100 ||^. The parameter A controls the flow of heat to and from the QC. 

The procedure for simulating a QC in the presence of dissipation is very similar to the one 
used in the pure quantum case: The only difference is that we solve the master equation (32) 
instead of TDSE (2). For A = the two completely different numerical methods used to 
solve these two equations (see |^ and [^] respectively) yield identical results. 



We have studied the effect of dissipation by performing calculations for the same cases 



as described in Section [IV A| . Our simple-minded model of dissipation is found to cause 
only decoherence. We don't find qualitative differences between the time-evolutions corre- 
sponding to the two different initial conditions. The presence of dissipation does not seem 
to affect the QPP. For A = 10~^ the value of the qubits in the final states is approximately 
1/2, i.e. all useful information is lost (see Fig. 0). Apparently our choice does not lead to a 
relaxation during the XX pulses. Although the dissipation processes incorporated through 
the use of master equation (32) lead to the correct thermal equilibrium state, there is no 
unique prescription to choose the bath Hamiltonian and the coupling between the bath and 
the spin system, i.e. the form of C. We are currently studying more realistic mechanisms 
for dissipation and will report on the results elsewhere. 



VI. FOUR-QUBIT QC'S 

In general practical applications of QC's will require complicated QA's. Therefore it is 
of interest to study how the instabilities due to imperfections in the operations affect the 
performance in such complicated QA's. In this section we use a copy (or swap) procedure 
in a four-qubit QC as one of the simplest examples. 

As discussed above, an implementation of Grover's search algorithm on a QC consist of 
two separate parts: 1) The preparation of the uniform superposition state \u) and encoding 
of the database information and 2) the search (query) of the database. The algorithms 
described above use the same two qubits for these two different tasks. First they are used 
to store the information contained in the database and second they are used to carry out 
Grover's algorithm to query the very same two qubits. Although this is sufficient to demon- 
strate the realization of 2-qubit QC hardware, from the point of view of computation this 
demonstration is of little use. Indeed, instead of organizing the computation in two steps 
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(bringing the two qubits into a state that reflects the content of the database and then using 
the same two qubits to perform the query) we could have written the information from the 
database directly into the qubits in a form which reveals the position of the item in a trivial 
manner. 

A computationally non-trivial demonstration of a QC running Grover's database search 
algorithm requires the database and quantum processor to have their own qubits. This 
obviously means using at least four qubits, two to hold the database information (in su- 
perposition state) and two to process the query. Therefore we will need an operation to 
copy the state of the database into the quantum processor. A schematic picture of how the 
calculation is organized is shown in Fig. |^. Let us now study how instabilities in the physical 
operations affect the outcome of this more complicated computational procedure. 

The copy operation transfers the amplitudes of an arbitrary linear combination of spin-up 
and spin-down of e.g. spin 1 to e.g. spin 2, initially in a state of spin up. More specifically 

Ci,2(a|0)i + 6|l)i)|0)2 = |0)i(a|0)2 + &|1)2). (34) 

In principle this can be done by a network of CNOT gates 0. One possible realization of a 
CNOT gate is the pulse sequence 

Ci,2 = riXiriXiFiX2F2/(7r)F2, (35) 
up to irrelevant phase factors, but there are many others. A fundamentally different sequence 



where 



/'(2a) 



i is 






Cl,2 - 


/ |oo) \ 


/ 


|10) 




|01) 




V|ii)^ 


V 



Y2X2Y2l'{7r) 



(36) 





cos a 
i sin a 






i sin a 
cos a 








( |oo) \ 







|10) 







|01) 


l) 







(37) 



A simple calculation shows that the two-qubit operation /'(2a) corresponds to the time- 
evolution of a two-spin XY model: 



I\2tJ, 



xy) 



(38) 



The other pulses in sequence (36) serve to remove unwanted phase factors. Note that (34) 
and (38) destroy the state of the database, a manifestation of the "no-cloning theorem" ||62|| . 
As Fig. ^ shows this is not really a problem as the state of the database after the copy 
operation took place is the same state that was used to initialize the state of the database. 

As before our aim is to address the fidelity of the results obtained by running a quantum 
algorithm on a hardware realization of a QC. Our QCE can simulate various candidate 
technologies and architectures (ways of interconnecting different units). However, to analyse 
the QPP it is at present sufficient to consider marginal extensions of two-qubit QC's. With 
this in mind we have chosen to "implement" a 4-qubit QC as two identical Ising models 



T QZ QZ 



hz{giSl + g2Sl 



(39) 
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In our numerical work we take for the values of the model parameters those corresponding 
to the chloroform molecule and we use selective resonant pulses to address the individual 
spins (see above). Thereby we assume, again for simplicity of analysis, that we can address 
pairs (1,2) and (3,4) separately (i.e. a pulse tuned to spin 1 does not affect spin 3 etc.). In 
the case where we use CNOT's to perform the copy-qubit operation we add to (39) the Ising 
interactions 

Hi,3 = -J.SISI (40a) 
H2A = -JzS^Sl (40b) 

or, in the case the XY-model is used to copy the qubits from the database into the query 
processor, we add to (39) 

Hl^ = -J,y{SlSl + SlSl), (41a) 
H',^, = ~J.y{S^2Sl + SlSl). (41b) 

It is obvious that these operations cannot be realized in terms of NMR technology. However 
this extension from a 2-qubit to a 4-qubit QC is sufficient to study in detail the stability of 
QA's. 



A. Resonant pulses 

In our numerical simulations on the QCE we take the same model parameters as those 
of the chloroform molecule (see above) for both the database and query machine. We use 
resonant pulses to address the single qubits. Thereby we assume that qubits 1 and 2 can be 
shielded from qubits 3 and 4 during the application of these pulses. 



The results of running Grover's search algorithm on the QCE are summarized in Table . 
In these calculations we used the pulses optimized in the sense discussed above (i.e. when 
used to run on a 2-qubit QC they yield the correct answers, independent of the initialization 
sequence). Nevertheless we observe that the implementation that uses the CNOT gates 
performs rather poorly, at the border of giving the wrong answers. As in the 2-qubit case, 
the fact that the duration of the pulses is no longer commensurate with the Larmor periods 
of the spins is one source of errors. Another mechanism that causes errors to accumulate is 
discussed below. 



B. Rotating field 



Finally, in an attempt to further reduce the uncertainty on the outcome of the quantum 
computations, we have repeated the simulations using resonant pulses of rotating fields 
More specifically, to apply a pulse to rotate for example spin 1, we add to the 



Hamiltonian a term of the form 



Hi{t) = hS] sin /it + hiS^ cos f^t, (42) 

where the frequency /i is tuned to the resonance frequency of spin 1 and the intensity of 
the pulse hi is determined by the desired angle of rotation. A simple calculation shows that 



14 



such a pulse transforms the single qubit in exactly the same manner as the corresponding, 
ideal transformation [^,|^. It is evident that this kind of optimization should improve the 
stability of the QA with respect to logically-allowed interchanges of operations. Further 
improvements can be made by reordering some of the elementary operations (recall that 
the sequence implementing e.g. a CNOT is not unique). The results of the simulations are 
presented in Table 0. Compared to the results of Table |^ there is one major change: The 
CNOT based implementation performs much better. 

In this implementation there are two sources of errors: The effect of the pulses on 
non- resonant spins and the interaction-controlled phase shift. In the case of the chloroform 
molecule there is a large difference in time scale (roughly a factor 10®) between the precession 
frequencies and the spin-spin interaction. Therefore, during the execution of an interaction- 
controlled phase shift, the spins make a large number of rotations about the z-axis. We 
have seen that it really matters for the stability of the calculation whether the angles of 
rotation are a multiple of 27r or not. However, in view of the difference in time scale, this 
implies that the duration of an interaction-controlled phase shift must be specified with a 
sufficiently high accuracy in order to recover the correct final results. 

VII. SUMMARY 

We have studied the difficulties that are encountered in programming quantum computer 
hardware. Taking Grover's database quantum search algorithm as an example, we ran vari- 
ous, logically identical versions of the algorithm on different, physically realizable quantum 
computer hardware. We demonstrated that the choice of the physical processes used for 
quantum computation has direct consequences in terms of programming the machine. 

We have shown that in the absence of dissipation non-ideal, physically realizable quantum 
computers operate in a regime of extreme sensitivity. This high sensitivity reflects itself in 
problems of programming quantum computers such that they perform correct calculations. 
These problems cannot be solved by error correction but may be reduced, and in some cases 
almost eliminated, at the cost of extensive, machine and program-specific optimization. 
A possible route to reduce this sensitivity may be to introduce some form of dissipation 
in a well-controlled manner. Dissipation can extend the region of stable operation of the 
quantum computer but also limits the time interval for quantum computer operation due to 
decoherence. We are currently exploring ideas along this line. 
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APPENDIX A: QUANTUM COMPUTER EMULATOR (QCE) 

In this appendix, we show an example of QCE. In Fig. | we show a snapshot of a window 
of QCE with a set of operations used to produce the results of Table | for the optimized 
pulses. 
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The left panel in Fig. ^ lists the fundamental operations, where MI stands for micro 
instruction. For example MI '-X1' means the operation Xi, etc. Mi's are defined by a 
table with interaction parameters, duration time, etc. This table appears when we double- 



click an MI. More details can be found in 37 . The set of Mi's for the case of resonant 



optimized pulses is stored under the name 'nmr-insta'. A QP represents a subroutine which 
consist of sequence of Mi's. For example QP 'f2' means the sequence of operations -F2, see 
Eq. (11). The MI 'initialize' simply prepares the state |00), MI 'inv-means' performs the 
inversion about the mean (the essential step in Grover's algorithm, see Eq. (13)), MI 'tau' 
corresponds to the operation /(vr) (see Eq. (8)), QP 'prepare' denotes the preparation of \u) 
while QP 'prapera' denotes the preparation of \u'), etc. The Mi's 'tau2', 'tauJnsta', and 
'wl' are not used in this example. 

The eight small windows show the programs corresponding to Uq, ■ ■ ■ U3 using \u) (upper 
level) and those using \u') (lower level), which are named QP 'gO' to QP 'g3prapl'. The 
result of running each program is shown in the grid below the QP. The expectation value of 
(1 — 5*^ )/2, 2 = 1 or 2, and a = x,y or z, can be read off from a particular cell by moving the 
mouse over the cell (not shown) . In Fig. ^ the approximate value can be determined from 
the darkness of the cell: Light gray corresponds to zero, dark gray to one. For example, the 
value of the cell {z, 1) in QP 'gO' is 0.027 and that of {z, 2) is 0.152, as given in Table | (Q; 
and Q2 respectively). For other cases, e.g. the non-optimized pulses, we use another set of 
Mi's ('nmr') with parameters given in . 
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TABLES 



TABLE I. Final state of the qubits after running the Grover's database search algorithm on an 

ideal QC {Qi,Q2), on a NMR-like QC {Qi,Q2), and on the same QC {Qi,Q2) using a different, 
but logically identical, initialization sequence = |n') instead of \u) = \u)). The optimized 
results, marked with a prime, have been obtained by changing the duration and intensity of the 
pulses tuned to the resonance frequency of spin 1 by a factor of two and 1/2 respectively. 



Uo\u) 




Ui\u) 


U2\u) 






Qi 




1 







1 


Q2 







1 




1 


Qi 0.028 




0.966 


0.037 




0.995 


Q2 0.163 




0.171 


0.836 




0.830 


Qi 0.955 




0.041 


0.971 




0.027 


Q2 0.031 




0.026 


0.971 




0.972 


Q[ 0.027 




0.972 


0.037 




0.964 


Q'2 0.152 




0.180 


0.847 




0.820 


Q[ 0.030 




0.969 


0.034 




0.965 


Q'2 0.022 




0.035 


0.977 




0.965 


TABLE XL Final state of the qubits after running the query part G of Grover's 


database search 


algorithm on an ideal QC {Qi,Q2) and 


on a NMR-like QC iQ'i,Q'2)- In all cases 


the input state 


|^'2) corresponds exactly to the case where the item 


is located in position 2 in the database. The 


queries G and G are logically identical to G. The differences in the outputs of the NMR-like QC 


are due to the internal quantum dynamics of the physical qubits used. 






G\^2) 


G\^2) 


G|*2) 


G^\^2) 


G^\^2) 


G^^2) 


Qi 

















Q2 1 


1 


1 


1 


1 


1 


Q[ 0.257 


0.154 


0.487 


0.908 


0.766 


0.578 


0.944 


0.944 


0.988 


0.938 


0.857 


0.995 
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TABLE III. Final state of the QC after running the Grover's database search algorithm on 
an ideal QC (QijQ-z) and an QC (Qi^Q-z) using hard, non-selective pulses. {Qi,Q2)- QC using 
hard non-selective pulses, and a different but logically identical initialization sequence (\u) = \u') 
instead of \u) = \u)). 



Uo\u) 




Ui\u) 


U2\u) 




U3\u) 


Qi 
Q2 




1 





1 




1 
1 


Qi 0.000 
Q2 0.000 




0.999 

0.005 


0.001 

0.999 




0.999 

0.996 


Qi 0.001 
Q2 0.001 




0.999 
0.002 


0.002 
0.999 




0.999 
0.998 


TABLE IV. Final state of the QC after running the Grover's database search algorithm on a 
4-qubit QC. Q^,Qa- ideal QC; Q^,Qa: NMR-like QC, using an XY-model time-evolution to copy 
the state of the database to the two-qubit QC that performs the query; Qs, Qa: NMR-like QC, 
using a sequence of CNOT to perform the copy operations. The values of the two database qubits 
are not shown. 


Uo\u) 




Ui\u) 


U2\U) 




Us\u) 


Qs 
Q4 




1 





1 




1 
1 


Q3 0.172 
Q4 0.135 




0.809 
0.187 


0.190 
0.865 




0.831 
0.814 


Q3 0.412 
Q4 0.113 




0.565 
0.218 


0.433 
0.888 




0.591 
0.783 


TABLE V. Final state of the QC after running the Grover's database search algorithm on a 
4-qubit QC, using pulses optimized to reduce the effect of imperfections of the operations and 
phase errors. Qa, Q4: ideal QC; QsjQa: NMR-like QC, using an XY-model time-evolution to copy 
the state of the database to the two-qubit QC that performs the query; Q3, Q4: NMR-like QC, 
using a sequence of CNOT to perform the copy operations. The values of the two database qubits 
are not shown. 


Uo\u) 




Ui\u) 


U2\U) 






Qs 
Q4 
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1 
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Q3 0.195 
Q4 0.195 




0.806 
0.194 


0.193 
0.805 




0.806 
0.806 


Qa 0.265 
Q4 0.197 




0.737 
0.193 


0.262 
0.804 




0.737 
0.807 
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FIGURES 




FIG. 1. Region of stable (black) and unstable (gray) operation of Grover's search algorithm 
when executed on a 2-qubit NMR-like QC. The x-coordinate (x = (^'l^fo)) gives the projection 
of the random input state 1^*') on the exact input state ^'o, i.e. for the case where the item is at 
position 0. The y-coordinate {y = |(00|G^'')|) discriminates between correct and false output. A 
confidence level of c = 0.7 on the input-output pair was used to determine if G yields the correct 
result. The horizontal (vertical) arrows mark the values of y (x) at which the input (output) state 
is equal to its ideal value. Note the absence of points in the regions near (—1, 1) and (1, 1). 
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FIG. 2. Time evolution of the qubits Q[ (solid line) and Q2 (dashed line) obtained by executing 
U2\u') on the physical model of an NMR-like QC, in the presence of dissipation (A = 10^^). The 
time intervals for each operation have been rescaled to make them look equal. 



22 




FIG. 3. Structure of a quantum program for Grover's search algorithm, using a 2-qubit QC 
to hold the database information and another 2-qubit QC to perform the search. 
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FIG. 4. Snapshot of the Quantum Computer Emulator showing a window with a set of op- 
erations implementing a quantum computer using resonant pulses and windows with quantum 
programs implementing Grover's database search. 
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